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ABSTRACT 

The form of depleted sulphur in dense clouds is still unknown. Until now, only two molecules, 
OCS and SO2, have been detected in interstellar ices but cannot account for the elemental 
abundance of sulphur observed in diffuse medium. Chemical models suggest that solid H2S 
is the main form of sulphur in denser sources but observational constraints exist that infirm 
this hypothesis. We have used the Nautilus gas-grain code in which new chemical reactions 
have been added, based on recent experiments of H2S ice irradiation with UV photons and 
high energy protons. In particular, we included the new species S„, H2S,, and C2S. We found 
that at the low temperature observed in dense clouds, i.e. 10 K, these new molecules are not 
efficiently produced and our modifications of the network do not change the previous pre- 
dictions. At slightly higher temperature, 20 K in less dense clouds or in the proximity of 
protostars, H2S abundance on the surfaces is strongly decreased in favor of the polysulfanes 
H2S3. Such a result can also be obtained if the diffusion barriers on the grains are less im- 
portant. In the context of the life cycle of interstellar clouds and the mixing between diffuse 
and denser parts of the clouds, the depletion of sulphur in the form of polysulfanes or other 
sulphur polymers, may have occurred in regions where the temperature is slightly higher than 
the cold inner parts of the clouds. 
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1 INTRODUCTION 



Contrary to other elements (see |jenkin^l2009h . the gas-phase abun- 
dance of atomic sulphur, in the diffuse medium, is observed to be 
constant with cloud density, around i ts cosmic value of ~ 1 .5 x 10"^ 
(compared to H) ( Sofia et al.l 1 19941) . In dense clouds (with den- 
sities above 10'' cm"'), the sum of the abundances of molecules 
containing sulphur is a small fraction of the cosmic abund ance of 
S (about 10 \ see lOhishi et"Zll 19921 : [Pickens et alj|200(3h . To re- 
produce such low S-bearing abundances, one needs to use a very 
small ini tial abundance of atomic sulphur of a few 10"* (see for 
example I Wakelam & Herbslll2008l) . This suggests that sulphur de- 
pletion occurs very efficiently and rapidly between the diffuse and 
dense phases of the cloud life leading to an unknown, still un- 
observed, reservoir of sulphur on grains. When an atom of sul- 
phur sticks to a grai n, models predicts that it would be hydro- 
genated to form H2S l lGarTod et alj|2007h . H2S has however never 
been detected in interstellar ices and its a bundance has been es- 
timated to be smaller tha n 5 x 10"* * (/H) (^van Di shoeck & Blakd 
Ll998; Jimenez-Escobar & Muifoz Ca ro 2011). Soli d DCS has been 
observed with an abundance of about 5 x 10"* jPalumbo et al.l 
[l993). SO2 may also be pres ent with a simila r abund ance according 



to 



Boogert et al. I ( Il997h and lZasowskietalJ ( l2009t) . Solid H2SO4 



has been proposed as a possible reservoir but no detection have 



confirmed this hypothesis jScappini et Z||2003|) . Fina lly, iron sul- 
fide Fe S has been observed in protoplanetary disks bv lKeller et al.l 
12003). The depletion of sulphur in diffuse clouds does however not 
follow the one of iron so that the iron left in dense clouds cannot 
account for such a depletion of sulphur. 

Laboratory experiments have r ecently brough t some new 
insights into the sulphur prob l em. [ Garozzo et al.l j2010l) and 
' Jimenez-Escobar & Mufioz Card ^201 A have independently re- 
cently published experimental results of the irradiation of H2S in- 
terstellar analogs ices by high energy protons and UV photons re- 
spectively. Both studies found that solid H2S was easily destroyed 
in favor of other species such as OCS, SO2, CS2, H2S2 etc. Based 
on these results, we have revisited the sulphur depletion problem 
in dense clouds introducing new mechanisms for the formation of 
polysulfanes (H2S,,), CS2 and sulphur polymers (S„) on the grains. 
The model, parameters and modifications of the chemical network 
are described in section |2] The results of our model are presented 
in section[3] We conclude in the last section. 



2 MODEL DESCRIPTION 

We have used the chemical model Nautilus, which computes the 
abundances of chemical species as a function of time in the gas- 
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Table 1. Formation reactions for tlie main S-bearing species 
Grain surface reactions 
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Table 2. S-polymers formation reactions on the grains (wliere tlie larger 
polymer formed is Sg)- 

Grain surface reactions 



JS;, + JS — > JS,j+i 
JS,j ^" JS,, ^ JS2f; 



phase and at the surface of interstellar grains. Starting from an 
initial composition, the model solves kinetic equations that de- 
scribe how species abundances are modified by the various reac- 
tions destroying or forming them. These reactions are typical gas- 
phase reactions (see lWakelam et al.l2010l . for a description of these 
processes), interactions between the gas and dust surfaces (ad- 
sorption of gas-phase species on grains and evaporation processes 
through various mechanisms) and surface reactions. Surface reac- 
tions stand for the reactions between two species on grain surfaces 
(this process implies that species move on the surface to meet), as 
well as dissociation of species on surfaces by external UV pho- 
tons and UV photons produced by the interaction with cosmic-rays 
(see section |23] |. Details on these pro cesses and on the eq uations 
used to describe them can be found in lHersant et al.l ( l2009h and in 
ISemenovetai] ( 120101) . 

2.1 Chemical network 

The chemica l network used for t his work is the same as the 
one used by iHincelin et alj ( l201lh . It contains 6142 reactions, 
of which 4394 are pure gas-phase reactions and 1748 are grain- 
surface and gas-grain interaction reactions. The model follows the 
chemistry of 458 gas-phase species and 195 species on grains. 
The gas-phase network is based on the osu-08-2009 database 
( |http://www.physics.ohio-state.edu/~eric/research.html| l with up- 
dates according to the recommendations from the experts of the 
KIDA database (http://kida.obs.u-bordeauxl.fr). 

To this network, we added reactions for a new species: CS2 
(as well as CSJ, HCS2 and CS2 on the surface). Since CO2 and 
CS2 are isoelectronic, we have based the new reactions for CS2 on 
those of CO2. Rate constants for these reactions have been taken 
as mean values among the rates of the same type of reactions. We 
also added new reactions of formation on the grains for S-bearing 
molecules for OCS, SO2 and H2S2. Formation of sulphur polymers 
S„ (n = 2-8) and polysulfane H2S3 on the grains was also added, for 
these molecules could be good candidates to a sulphur refractory 
residue (Steudel 2003). This extends the network to 666 species 
(461 for the gas-phase and 203 on the grains) and 6215 reactions. 
The reactions that we added to the network are listed in Tables [T] 
[2l|3]and|4] and the full network is available on the KIDA website 
(http://kida.obs.u-bordeau xl.fr/modelsl l. 



[h] 

Table 3. Formation and destruction reactions of the polysulfanes on the 
grains (where the larger considered polysulfane is H2S3). 

Grain surface reactions 



JH2S2 + JS ^ JH2S3 
JHS2 + JHS — ^ JH2S3 

Grain surface CR photodissociation 



JH2S3 JHS + JHS2 

Grain surface UV photodissociation 

JH2S3 — » JHS + JHS2 

2.2 Model parameters 

For the simulations of the chemistry in dense clouds, the gas and 
dust temperatures are set to 10 K. The H density is 2x lO"* cm"\ the 
visual extinction is 30 and the total H2 cosmic ray ionization rate is 

1.3 X 10"" s"'. All elements are initially in the atomic form and in 
the gas-phase (assuming that dense clouds are formed from diffuse 
medium), except for hydrogen, which is supposed to be entirely in 
the form of H2 in the gas. Initial abundances are listed in Table[5] 
Elem ental abundance o f sulphur is at its cosmic abundance of 1.5 x 
10-5 ( ISofiaetal.lll994h . 

2.3 Model parameters for surface reactions 

For the formation of molecules on the surfaces, there are two 
key processes for which we give some details: the dissociation of 
species at the surface of the grains by photons induced by cosmic- 
ray particles and the reactions occurring at the surfaces. In our 
model, the efficiency of these processes depends on a number of pa- 
rameters that we describe now. The direct interaction of cosmic-ray 
particles with chemical species produc es mostly ionization of H 2 
and atoms (.Cravens & Dalg arno 197^: IPrasad & Huntres3ll980l) . 
Ionization of H2 by cosmic-rays produces high energy electrons 
that excites H2. By de-excit ation, H2 th en produ ces a UV radiation 
field inside dense clouds ( jPrasad ^Tarafdar 1 9831 ; lOredel et al.l 
Il989l) . The dissociation of molecules by these photons are usually 
called cosmic-rays induced photodissociations. Such processes are 
also supposed to be efficient on the molecules contained in grain 
ices and their efficien cies are t aken to be the same as for gas- 
phase molecules as in lRuffle"& Herbst (2000). Cations produced 
by photo-ionizations are assumed to recombine on the surfaces of 
the grains, which are mainly negatively charged. Rate coefficients 
for these reactions are computed by the formula: kcjtp = acRp( 
(in s"'). ^ is the total H2 cosmic-ray ionization rate and ocrp 
specific for each species (see Gredel et al. 1989). Table |6] lists 
the acRp used in our model for the main S-bearing molecules, 
which are from the OSU databasfl For H2S3 on the surfaces, we 
have used the same photodissociation rate as for H2S2. The prod- 
ucts of the photodissociation of H2S2 on the surfaces have been 
ch anged compared to the gas-phase based on the new experiments 
bv I Jimenez-Escobar & Mufioz Card ( l201lh . 

Reactions on grain surfaces between two chemical species are 
computed using equation 5 from lHasegawa et "all d 19921) : 

Rij = KijiRdiffj + Rdiffj)NiNjnd (1) 

with Rj j the rate (in cm"^ s"') of reaction between species ; and 

' http://www.physics.ohio-state.edu/~eric/research.html 
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Table 4. Gas-grain CS2 chemistry 

Grain surface reactions 



Table 5. Initial abundances fsee lHincelin et alj201lL for discussion on these 
values) 
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Gas phase UV photodissociation 
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Element 


Abundance (/H) 


H2 


0.5 


He 


9 X 10"- ^ 
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Mg+ 
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" IWakelam & HerbsllfcOOSi) 

^ IjenkinslfcOOgI) 

" lHincelinetai]j201lh 

See text 

■= Low metal elemental abundances, (see lGraedel et al.l 19821) 



Table 6. acRp for the main S-bearing species on the grains. 



CR induced photodissociation reaction 
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i at the surface of the grains, /c,j is the probability for the reac- 
tion to occur (one for exothermic reactions). A', and Nj the num- 
ber of species i and j at the surface, is the number density of 
grains (in cm"^). Rdiff are the diffusion rates of t he species at the 
surface of the grains (see equations 2 to 4 from iHasegawa et alj 
[T992). This diffusion rate is proportional to exp(-Ei,/kTj) with Tj 
the dust temperature and £4 the potential energy barrier re quired to 
move from one site to the other on the surface of the grains jWatsonl 
Il976l) . This last parameter is uncertain and is usually taken to be 
a fraction of the binding energies (Ed) of the species to the sur- 
face. Various values for this fraction have been used over the years, 
typically ranging from 0. 3 jWatsorJI 19761: iTielens & Allamandolal 



ll987l : lHasegawa et alj[l992l) to 0.77 dRuffle & Herbstll2000l) . Here 
we use a mean value of 0.5 as in iGarrod & Herbslltoel) . 



3 MODELLING RESULTS 

3.1 Under dense cloud conditions 

The model was used with the parameters described in the previ- 
ous section. The abundances (compared to total H density) for the 
main S-bearing molecules in the gas-phase and on the grains are 
shown in Figs.[T]and|2]as a function of time. In the gas-phase, the 
predicted abundances of CS, SO and H2S are much larger than the 
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H,S 



10^ yr 



102 103 lO-* 105 106 10^ 108 
log(time [yr]) 

Figure 1. Gas-phase abundances of the S-bearing species relative to H as a 
function of time for a gas and dust temperature of lOK. 
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Figure 2. Grain surface abundances of the S-bearing species relative to H 
as a function of time for a gas and dust temperature of 10 K. Red horizontal 
dash ed line represents the upper li mit on H2S abundance observed in grain 
ices ivan Dishoeck & Blakdll9 



ones observed in the dense cloud L134N (North Peak). A compar- 
ison between the predicted (at 10* yr) and observed abundances is 
listed in Table |7] Observed abundances from several authors are 
given to see the level of uncertainties in the observations. In paral- 
lel, sulphur, which is initially in the atomic form in the gas-phase, 
depletes efficiently on grains and is hydrogenated to form H2S. Hy- 
drogen sulfide on grains is then the main reservoir of sulphur and 
its abundance predicted by the model is several orders of magnitude 
larger than the observational upper limit, similarly to previous mod- 
elings (see for instance Garrod et al. 2007). The other S-bearing 
molecules efficiently produced on the grains are H2S2 and CS2 but 
with abundances of 1.1 x 10"^ and 2.2 x 10"'" respectively at 10* yr. 
At low temperatures as 10 K, the only efficient grain surface reac- 
tions are hydrogenations, which reform H2S when dissociated. The 
modifications that we have done on the chemical network do not af- 
fect the results here, except that OCS is less produced on the surface 
in favor of the newly introduced species CS2. For other S-bearing 
molecules to form on the grains, the diffusivity of the species on 
the grains have to be increased. 
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Figure 3. Grain surface abundances of the H2S, H2S2, H2S3 and Sg as a 
function of temperature at two dilferent times (10* yr on the top and lO' yr 
at the bottom). 



3.2 Sensitivity to the model parameters 

We now test sensitivity of the model results to some of the key 
parameters. 



3.2.1 The gas and dust temperature 

Figure |3] shows the grain surface abundances of the H2S, H2S2, 
H2S3 and Sg at two different times (10* and 10' yr) for tempera- 
tures (gas and dust) from 10 to 70 K. At 20 K, half of the sulphur 
goes into H2S3 on the surface after 10* yr and the abundance of 
H2S is below 10 ' after 10' yr (see also Fig.©. In this model, CS2 
is produced efficiently and its abundance is similar to that of H2S 
for times between 10* and 10' yr. At a temperature of 30 K, the 
abundance of H2S on the grains is at maximum 10"^ (at 10* yr). Its 
gas-phase abundance is always smaller than 10"'". Above 30 K, the 
abundance of S-bearing molecules on the grains decreases because 
the evaporation temperature of atomic sulphur is reached. At these 
low densities (lO* cm"'), the main form of sulphur is then neutral 
atomic sulphur in the gas-phase. The molecule Sg is efficiently pro- 
duced by temperatures between 30 and 50 K. 



3.2.2 Energy barrier between two grain adjacent sites 

Another parameter that inffuences directly the diffusion of the 
species on grains is the energy barrier between two grain adjacent 
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Table 7. Abundances relative to H observed in L134N (North Peak). 
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Element Abundances (/H) of the main gas-phase S-bearing species. 
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Swade(1989) 


Ohishietal. (1992) 


Minh et al. (1989) 


Dickens et al. (2000) 
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4 X 10"** 


3.4 X 10"'" 


5x 10"'° 




8.3 X 10"'° 


so 


1 X lO-** 


0.8-2.1 X 10"'° 


1 X 10"** 




1.6 X 10"' 


S02 


1 X 10"'° 


< 1.4 X 10"'° 


2 X 10"' 




< 8x 10"'° 


ocs 


1 X 10"' 




1 X 10"' 






H2S 


4 X 10"* 




4x 10"'° 


1.5 X 10"' 





" from two dilferent molecular lines. 




102 103 104 105 106 10^ 108 

log(time [yr]) 

Figure 4. Grain surface abundances of the S-bearing species relative to 
H as a function of time for a 20K gas-grain model. Red dashed hor- 
izontal line represents the u pper limit on H2S observed in grain ices 
ivan Dishoeck & Blakdl 19981) . 
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Figure 5. Grain surface abundances of the S-bearing species relative to H 
as a function of time for Et = 0.3Eo (other parameters are the same as for 
Fig.0. 

sites £■/,. As mentioned in the section 2.4, the value of this energy 
barrier is far from being well known. To test its importance, we 
have run the model described in section 2.3 (dust and gas temper- 
ature are 10 K) with an Eb of 0.3£d (see section 2.4). The model 
results are shown in Fig. |5] The destruction of H2S on the surfaces 
is very efficient in that case in favor of the production of CS2 and 
H2S3. After a few 10^ yr, H2S3 and H2S2 are the main S-bearing 
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Figure 6. Grain surface abundances of the H2S, H2S2, H2S3 and Sg as a 
function of temperature at two dilferent times (10* yr thin line and lO' yr 
thick line) for a H density of 2 X lO' cm"^^. 



species on grains and the prediction abundance of H2S is below its 
observational limit. In that case, CS2 is not produced efficiently. In 
all the models shown in this section, the formation of pure sulphur 
polymers S„ is very low and the predicted abundances are below 
the minimum values shown on the figures. 



3.2.3 Other parameters 

At a larger density of 2 x 10^ cm"^, the results at 10 K are very 
similar to the lower density one. Between 20 and 60 K, solid H2S3 
remains with a very high abundance similar to the neutral atomic 
sulphur in the gas. The molecule Sg is efficiently produced for 
temperatures between 30 and 60 K. At 70 K, gas-phase S is the 
reservoir of sulphur. A figure similar to Fig. [3] for a density of 
2 X 10'' cm"^^ is given in Fig.|6] 

The results presented here are not strongly sensitive to the cosmic- 
ray ionization rate. A ten times larger L,, with a temperature of 
20 K, accelerates the destruction of H2S on the surface and favors 
the production of H2S2. H2S3 remains however the reservoir of 
sulphur after 5 x 10' yr. A ten times smaller ^ only produces a 
smaller abundance of H2S2, below 10"'° for times between 10* 
and 10' yr. 



4 CONCLUSIONS 

Chemical models that take into account the formation of molecules 
on interstellar grains predict that sulphur would stick on the 
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grains at high density and would be hydrogenated to form 
H2S. The non detection of this species in interstellar ices tends 
to indicate that solid H2S cannot be the reservoir of sul- 
phur. The recent experiments f rom i Garozzo et alj jZOlOl) and 
I Jimenez-Escobar & Munoz Carol bOllb showed that H2S is eas- 
ily dissociated by high energy particles and UV photons, and that 
other molecules such as CS2, polysulfanes (H2S,,) and S„ would be 
formed on the surfaces. Using a gas-grain model, in which new 
reactions have been introduced for CS2, H2S,, and S„, we show 
that we are able to produce large quantities of solid H2S3 and very 
low abundance of solid H2S but only with temperatures higher than 
10 K, which are probably more representative of star forming re- 
gions. The key parameter here is the diffusivity of the radical HS 
on the grains, which is limited at very low temperature. Significant 
amount of Sg is produced at temperatures between 30 and 50-60 K 
for densities between 2 x lO'' and 2 x 10' cm"'. The results of our 
model, as for any chemical model, depends on the chemical re- 
actions considered. More experience on the diffusion of S-bearing 
species on surfaces at low temperature and the formation of these 
chains would be a precious ally. 

The possibility to observe products of the destruction of 
solid H2S,, or S„ in diffuse or warm medium could bring more 
stones to this problem. H2S2 is a symmetric molecule and the pre - 
dicted lines in the millimeter range are weak (lBehrendetalJl990h . 
Search for these line s, as well as S3 and S4 jThorwirth et alj 20051 : 
lOottlieb et al.l l2005'). in the s pectral surve y of the low mass pro- 
tostar IRAS 16293-2422 ( Caux et al.ll201 ih . was unsuccessful (E. 
Caux private communication). The results of this modeling can be 

i iut in the con text of the larger picture of interstellar dust cycle 
Tielenj[l998h and the fact that the depletion of sulphur into H2S,, 
may have happened at an early phase of the cloud history when the 
dust temperature was larger than what is observed at the centre of 
these cold cores. It is unlikely that simple diffusion mechanisms at 
the surface of the grains could form pure sulphur polymers (S„). 
The observation of OCS and SO2 in interstellar ices may indicate 
however that other mechanisms could take place or that surface re- 
actions are much more efficient than current models assume. 
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